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Abstract 

We present the study of the influence of the leading contribution 
from the kinetic energy of the Skyrme Lagrangian on the scattering 
of low energy nucleons. This classical kinetic energy for the 2-body 
system is computed using the product ansatz and is correct for low 
energy, well separated Skyrmions. We quantize the rotational degrees 
of freedom of the Skyrmions in order to give them the right spin and 
isospin values. We are then able to compute analytically the scattering 
angles for some angular momenta polarizations. 

1 Introduction 

In this work|j|] we present some of the implications of the Skyrme 
model on the scattering of two nucleons. We will first start with a 
short introduction to the Skyrme model and its simplest solution, the 
Skyrmion. Then in the next section we study the problem of the scat- 
tering of two Skyrmions using the product ansatz and present the con- 
straints imposed by this parametrization to the scattering processes. 
The expression of the energy of the two Skyrmion system is then writ- 
ten as an expansion in inverse powers of the separation between the 
particles. We then briefly describe the approximation methods used 
to find and solve the equations of motion of the system, namely the 
method of Lagrange and the method of variation of constants. In 



2 



order to study the scattering of nucleons, we have to quantize the re- 
maining degrees of freedom of the system, namely the rotation of each 
Skyrmion so as to give them the proper spin and isospin. This is done 
with the Bohr-Sommerfeld quantization technique. We are then able 
to compute the scattering of two nucleons for certain spin and isospin 
polarizations. 

2 The Skyrme model and the Skyrmion 

The Skyrme mode was introducedj^ in the 60's and is described by 
the following Lagrangian: 

^sk = -Y HU^d^UU^d^'U) + ^ tr{[U^d^U, U^d.Uf) (1) 

where U{x^) is an element of SU{2) which represents a pseudoscalar 
massless particle, namely the pion. and e are phenomenological 
parameters related to meson decay and low energy scattering, and are 
theoretically computable using QCD. 

The simplest stable solution of the Skyrme model is the Skyrmion: 

Us{x) = e^-^M-^ (2) 

where r = (r^,r^,r^) are the three Pauli matrices and the function 
F{r) is a 

monotonous function of r which decreases from vr at the origin, to 
at infinity. F{r) has to be computed numerically but falls off at 
large distance like n/r'^. Since the Skyrmion field attains the SU{2) 
identity when r ^ oo, we can consider the three dimensional space 
E? to be topologically compactified into a 3-sphere 5^. Us{r) then 
defines an application from an to the manifold of the SU{2) 
group. This type of application is classified into disjoint homotopy 
classes labeled by the winding number of the first 3-sphere onto the 
second, A^, defined by 

X 

247r^ J 

An application belonging to one class cannot be deformed continuously 
into another: it is topologically stable. The Skyrmion has winding 
number 1, rendering it stable against deformation into the vacuum 



defined by = 0. It was Skyrme who first interpreted as the 
baryon number. 

The Skyrmion is a static solution. In order to set it into motion, 
we shift its position by the time dependent vector R{t) and conjugate 
it by the time dependent SU{2) matrix A{t): 

Ui{r,t) = A{t)Us{x-R{t))A\t). (4) 

We then replace Ui (r , t) into the Skyrme Hamiltonian density and 
integrate over all space. This gives the following expression for the 
energy: 

Ei=M + yi?2 + 2Kr{A)C''{A) (5) 

where we defined the generators of the right and left action of the 
SU{2) group: 

7^"(^) = -1 tr[T''A^doA] (6) 

C'iA) = tT{T^d^AA\ (7) 

E\ has the familiar form of the energy of a translating and spinning 
spherically symmetric rigid body if M, of the order of 1 GeV, is in- 
terpreted as the mass, and A as the moment of inertia. The last term 
is the rotational energy, with L"'{^A)L°'{^A) equal to the square of the 
angular velocity. C^i^A) is usually interpreted as the tensorial part of 
the isospin of the Skyrmion, and 7^" (A) as its spin. 

It was shown by Atkins, Nappi and WittenQ, that by quantizing 
the rotational degrees of freedom A{t) the characteristics of the nu- 
cleon could be reproduced with an error varying between 10% and 
30%. This guides us towards investigating the predictions of the 
Skyrme model in the baryon number 2 sector. 



3 The Skyrmion-Skyrmion interaction 

We will now attempt to compute the classical Lagrangian describing a 
pair of Skyrmions in order to study their scattering. To achieve this we 
need a parametrization of the pair of Skyrmions. The simplest is the 
so-called "product ansatz" where the Skyrmion fields are multiplied 
together: 

U2{r,t) = Ui{Skyrinion 1) Ui{Skyrmion 2) 

= A{t)Us{x - Ri{t))A\t) B{t)Us{x - R2{t))B\t). ^ ^ 



A and B are SU{2) matrices representing the iso-orientation of the 
Skyrmions and i?i and R2 their position. This parametrization is not 
completely general and we have to restrict the initial conditions in 
order to obtain physically correct results. 

First, it is a well known result from numerical studies of the Skyrme 
model IQ] that Skyrmions deform when they come close to each other. 
In fact, the bound state of two Skyrmions, putatively the deuteron, 
has the form of a torus. These deformations cannot be described 
by Eq. ^, let alone the product ansatz, so we have to consider only 
configurations where the distance d = \Ri(t) — ^2(^)1 is large. 

Second, it was shown by N.S. Manton||5| that only low energy soli- 
ton systems can be described with a finite number of degrees of free- 
dom. We refer the interested reader to his article and will only say 
here that it is physically evident that the number of degrees of freedom 
of a system rises with its energy. Then if we keep the energy of the 
system sufficiently low, the system might excite only a finite number 
of modes. If we wish to describe our system with the 12 degrees of 
freedom of the product ansatz (6 to define the positions of the par- 
ticles, and another 6 for their iso-orientations), we have to consider 
only low energy scattering. 

We now give the expression of the energy for the pair of Skyrmions, 
in the center of mass frame, as an expansion in inverse powers of the 
relative distance d which is assumed to be large at all times: 

T = \Md^ + 2A(£'^(A) + £"(5) C'iB)) 

j^A^iac^jbd 7^c(^) n'^{B) - d'd^) Dab{A^B) + 0(l/(i2). 

(9) 

A = Ittk'^/^ is a result of the integration over space, d = d/d and 
DabiG) = 1/2 tr[r''Gr''Gt] is the 3 x 3 representation of the SU{2) 
matrix G. This term, induced by the kinetic energy of the Skyrme 
Lagrangian, was found by us[^ and independently by B.J. Schroers[0]. 

We note that this term has a structure similar to the spin-spin and 
tensorial terms found in the traditional nuclear potentials. Also, the 
leading contribution from the potential part of the Skyrme Lagrangian 
behaves as thus is neglected to leading order of our expansion of 
the energy. The kinetic energy defines a metric on the tangent space 
of the modulii space, the motion implied by the Lagrangian simply 
follows the geodesies of this metric. This is what we call the geodetic 
approximation. 



4 Equations of motion of the system 
and quantization 



We now have to find the equations of motion corresponding to the 
Lagrangian of Eq. ^. Because of the very comphcated form of the 
interaction term, we have to resort to various approximation meth- 
ods. We first choose the observables of the system, meaning, what 
we need to know from the scattering in order to possibly compare the 

results with experimental data. We choose the relative velocity d, the 
tensorial part of the spins TZ"'{A) and TZ°'{B), and isospins C"'{A) and 
£"(B). 

The approximation method of Lagrange |8| is perfectly suited to 
compute the time derivatives of the observables without having first 
to derive the equations of motion of ^, -B and d. The method uses the 
Poisson brackets formalism and is based on the following principle. Let 
us consider the degree of freedom g*. Conjugate to is the momentum 
which can be written as an expansion in powers of 1/d: 

p'=ph + Sp\l/d). (10) 

Then the Poisson brackets, which are functions of the and the p*, 
also form an expansion in powers of 1/d. If we denote by C'^ the 
observables of the system, the time derivative of is 

where we have divided the Hamiltonian into a free Hq (of order 1/d^) 
and interaction Hj part (of order 1/d and higher). By dividing also 
the Poisson bracket into a free and interaction part, we get 

±C'' = {C^Ho}i + {C\Hi}o (12) 
plus higher order terms. By computing the various Poisson brackets 



involved, we get the following set of 5 coupled equations: 
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fnf{A)Dab{A^B) + e''''fDfb{A^B)n^{A) 
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e^'^Snf{B)Dab{A^B) + e^'^f Daf{A^ B)1Z'^{B) 

(13) ■ 

where the dots represent other complicated terms. This is still too 
complicated to solve algebraically and we have to resort to further ap- 
proximations. We choose the approximation method of the variation 
of constants P]. Let us consider as an example the following equation: 



(14) 



and treat it not as an ordinary differential equation but more like an 
iteration equation: 

±xit)^fiMt),t) (15) 

where xo(^) is a trial function. Then if the trial function xq^I) is 
close to the true solution x{t) or if /(xo(i),i) is small, then a good 
estimation of the solution of the equation is: 



+ 



+ 



xit) 



f{xo{t'),t')dt'. 



(16) 



In our case, xo{t) is the free value of the C'^, so without the interac- 
tion term. The approximation method should give accurate results if 
the C'^ are slowly changing quantities. This is compatible with our 
working hypothesis. This method decouples the 5 equations of mo- 
tion and trivializes the integration over time of the observables: all we 

have to do is replace in each right hand side the free values for d, d, 
TZ''{A), TZ"-{B), CiA), CiB), A and B, and integrate over time each 
equation separately. 



We will not consider here the time evolution of the spins and 
isospins. We will only say that spin-flips and isospin-flips seem to oc- 
cur, which are consistent with spin 1 and isospin 1 particle exchange. 
We will only consider from now on the equation for the relative mo- 
mentum p = M/2d: 



_d i _ _A 



5'id''+y^d'+5^''d^-2,d^d^dl' 



e''"'€^^'^n\A)n'^{B)DaM^B) 



(17) 

in order to extract the scattering angles. 

To summarize, we now only have the equation for ^to solve. This 
is done by replacing in the right hand side of the equation the free 
values of d and the matrices A and B which deflne also the spin and 
isospin of the Skyrmions, and to integrate over time. 

So far we have only talked about classical scattering. If we want 
to study the physics of nucleons, we have to quantize the remaining 
degrees of freedom, i.e. the rotation of each Skyrmion. This means 
finding the matrices A and B which give the Skyrmions spin and 
isospin 1/2. The Bohr-Sommerfeld method is perfectly suited to this 
task. In the well known problem of the hydrogen atom, the sum of 
the action variables 

Ji= f Pidqi (18) 



is set equal to a multiple of the Plank constant h. This fixes the 

allowed values for the angular and linear velocities, thereby quantizing 
the energy and angular momenta of the states. In our case, the same 
program is applied to the time derivatives of the angles with which 
the matrices A and B are expressed. Here are the main steps of this 
procedure. 

Let us consider the Lagrangian for a single Skyrmion and express 
the matrix A of its iso-orientation as a function of 3 Euler angles a, 
/? and 7: 

A = ao + ia-T = e-^°^3/2 g-i/3T2/2 ^^g^ 

giving 

L = -M+^A [a^ + /^^ + 7^ + 2a7 cos P]. (20) 

The equations of motion are satisfied if a and 7 are linear functions of 
time t and if P is equal to or tt. The action variables for the angles 



a and 7 are then readily computed and give: 



J a = §Pada = 27rA[d + 7cos/3] = -A7rAC^{A) = -27r/3 

= fp^dj = 2^A[7 + d cos p] = -47r7^3(^) = 27rS3 ^ ^ 

using Eq. ^ and Eq. |^ and following the convention used by Adkins|^ 
et al. for the component of the spin and isospin along the z axis 



53 = 2A7^3(A) 
/3 = -2KC^{A) 



3.4^ (22) 



respectively. The Bohr-Sommerfeld quantization condition states that 
the sum + is equal to an integer multiple of h. This fixes the 
value of the angular speeds d and 7 so as to give, in our particular 
case, spin and isospin ±1/2. The expression for the matrix A giving 
the proper quantum numbers is then obtained by extracting the values 



of a{t) and 7(t) from Eq. 21, after choosing the value of fi which is 



either or vr. If /3 = 0, solving for the angles a and (i gives the matrix 

— gjM+(/'o)T^/2_ ^23) 

uj is the quantized quantity here, and has an absolute value of roughly 
100 MeV for the nucleon. This matrix represents the state \p [> if 
a; > and |n t> if a; < 0. li (5 = n, then we obtain 

which represents the state \p ]> if u > and [n J.> if w < 0. 

Before going further we note that the time dependence of the ma- 
trices A when /3 = or /3 = tt are very different. This will have 
important consequences on the nucleon-nucleon scattering. 



5 Nucleon-Nucleon Scattering 

We now have all the necessary tools to solve our problem. To solve 



Eq. |17| for a particular scattering process, we only replace the corre- 
sponding matrices A and B in the right hand side of the equation, as 
well as the free value of d: 

d={vt,-f,0) (25) 



which describes a relative particle traveling along the x axis but with 
an impact parameter 7 along the y axis. To respect the restrictions 
underlined in section 3, is chosen small to guaranty low energy, 
and 7 large, so as to keep the particles far from each other. In what 
follows, we choose z as the axis of polarization of the angular momenta. 
This gives 2-dimensional trajectories. Other polarization axes give 
complicated 3-dimensional motion best studied numerically. 

Since in our approximation scheme the spins and isospins are con- 



stant, the only time dependence on the right hand side of Eq. |17 
comes from d, d and Dah{-A'^B). It is then predictible that scattering 
processes separate into two cases, depending on whether or not A'^B 
is time independent. 

Dabi^L^ B) time independent: 



cos(2(5)/ (26) 
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[d^ + Ar'^f ■ d - Qd}'{r ■ df] 
r'^ = (— sin((5), cos(5), 0) 

[d^ + Ar^r ■ d-Qdf'if ■ df] 

= (sin((5), -cos(5),0) 



(27) 
(28) 



where 5 = (f)Q — cf)^ . The scattering angle defined as the angle between 
the initial {t = —00) and final {t = +00) momenta is the same for all 
these different processes and readily computed: 

cost' - |p-(+oo)||p(-^)| 

^ M-iv^ 1 (29) 

Dah{A^B) time dependent: 

d.pk = ^[dk+ 4r>'r ■ d - &d^ (f • df] , s 

^ ' ' f^ = {- sin(2u;t + e), cos(2a;t + e), 0) ' 



pin[ 



= (- sin(2u;t + e), - cos(2a;t + e),^0)'^ 



where e = 4>q + (j)^ . For those scattering states, the right hand side 
of the equations oscillates very quickly and the scattering angle is 
suppressed by the very small factor 



since u and 7 are large while v is small. 

We underline that this is the first analytical calculation of nucleon- 
nucleon scattering from essentially first principles, without recourse to 
ad hoc models or potentials. We emphasize that the Skyrme model is 
in principle derivable from QCD and /tt and e are, as such, calculable 
parameters and, in that sense this is also a QCD calculation. To 
calculate the classical scattering cross-section we need to compute the 
scattering for all different polarizations relative to the initial scattering 
plane. This would comprise a different project which would probably 
be best achieved by numerical methods. Therefore we are unable at 
this point to make a direct comparison with the experiment. 

A peculiar aspect of our treatment is the strong dependence of 
the scattering angle on the variables 5 and e which are just the phase 
lag between the rotation of the Skyrmions at some time t before the 
scattering. These act as hidden variables, and have to be measured 
after the scattering or averaged over. 
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